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ABSTRACT 

These last years ground-based astronomy has been looking towards Antarctica, espe- 
cially its summits and the internal continental plateau where the optical turbulence 
appears to be confined in a shallow layer close to the icy surface. Preliminary mea- 
surements have so far indic ated pretty good value for the see ing above 30-35 m: 0. 36" 
(|Agabi et al.ll2006[ ). 0.27" ([Lawrence et all I2004D and 0.3" (iTrinauet et all I2008D at 
Dome C. Site testing campaigns are however extremely expensive, instruments pro- 
vide only local measurements and atmospheric modeling might represent a step ahead 
towards the search and selection of astronomical sites thanks to the possibility to re- 
construct 3D maps over a surface of several kilometers. The Antarctic Plateau 
represents therefore an important benchmark t est to evaluate the p ossibility to dis- 
criminate sites on the same plateau. Our group (Hagelin et al. '2008) has proven that 
the analyses from the ECMWF global model do not describe with the required accu- 
racy the antarctic boundary and surface layers in the plateau. A better description 
could be obtained with a mesoscale meteorological model. The mesoscale model Meso- 
Nh has proven to be reliable in reproducing SD map s o f optical turbulence above mid- 
latitude as tronomical sites (Masciadri et alj Il999al lbl: iMascia dri and Ja bouilld l200ll : 
iMasciadrT et al. 2004). In this paper we study the ability of the Meso-Nh model in 
reconstructing the meteorological parameters as well as the optical turbulence above 
Dome C with different model configurations (mo no model and grid- nesting). We con- 
centrate our attention on the model abilities in reproducing the optical turbulence 
surface layer thickness {hsi) and the integral of the Cf^ in the free atmosphere and in 
the surface layer. It is worth to highlight that these are the first estimates ever done 
so far with a mesoscale model of the optical turbulence above the internal Antarctic 
Plateau. 
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1 INTRODUCTION 

The internal Antarctic Plateau is, at present, a site of po- 
tential great interest for astronomical applications. The ex- 
treme low temperatures, the dryness, the typical high alti- 
tude of the internal Antarctic Plateau (more than 2500 m), 
joint to the fact that the optical turbulence seems to be 
concentrated in a thin surface layer whose thickness is of 
the order of a few tens of meters do of this site a place in 
which, potentially, we could achieve astronomical observa- 
tions otherwise p ossible only by space. In spite of the excit- 
ing first result s (iLawrence et al.l l2004l : lAristidi et al.1 l2005l : 
iTrinquet et al. I l2008l ) the efi'ective gain that astronomers 
might achieve from ground-based astronomical observation 
from this location still suffers from serious uncertainties 
and doubts that have been pointed out in previous work 
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iHagelin et al.ll2008l : IStoesz et al.l l2008V A better estimate of 
the properties of the optical turbulence above the internal 
Antarctic Plateau can be achieved with both dedicated mea- 
surements done simultaneously with different instruments 
and simulations provided by atmospheric models. Simula- 
tions offer the advantage to provide volumetric maps of 
the optical turbulence (C|^) extended on the whole internal 
plateau and, ideally, to retrieve comparative estimates in a 
relative short time and homogeneous way on di fferent places 
of the plateau. In a previous paper (Hagelin et al.l 



: places 

iHooi) 



our group performed a detailed analysis of the meteorolog- 
ical parameters from which the optical turbulence depends 
on, provided by the General Circulation Model (GCM) of 
the European Center for Medium-range Weather Forecasts 
(ECMWF). In that work we quantified the accuracy of the 
ECMWF estimates of all the major meteorological param- 
eters and, at the same time, we pointed out which are the 
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Figure 1. Orography of Antarctica as seen by the Meson-Nh model (polar stereographic projection) in (a) the monomodel simulation, 
with horizontal resolution AX=100 km; (b) zoom of (a) above the Dome C area (the black square represents the same area as in (e)); (c) 
the largest domain of the grid-nested simulation with AX=25 km; (d) the second domain of the grid-nested simulation with AX=h km; 
(e) the innermost domain of the grid-nested simulation with /\X=1 km. The dot labeled 'C is located at the Concordia Station. The 
dot labeled 'A' is the Dome A. Altitude in meter (m). 



limitations of the GCMs. In contexts in which the GCMs 
fail, mesoscale models can supply more accurate informa- 
tion because they are conceived to reconstruct phenomena 
that develop at a too small spatial and temporal scale to 
be described by a GCM. In spite of the fact that mesoscale 
models can attain higher resolution than the GCMs, some 
parameters, such as the optical turbulence, are not explicitly 
resolved but are parameterized, i.e. the fluctuations of the 
microscopic physical quantities are expressed as a function 
of the corresponding macroscopic quantities averaged on a 
larger spatial scale (cell of the model). For classical mete- 
orological parameters the use of a mesoscale model should 
be useless if GCMs such as the one of the ECMWF could 
provide estimate with equivalent level of accuracy. For this 
reason the HageUn et al. paper (2008) has been a first step 
towards the exploitation of the mesoscale Meso-Nh model 
above the internal Antarctic Plateau. In that study we re- 
trieved an exhaustive characterization of all the meteoro- 
logical parameters from the ECMWF analyses (wind field, 
potential temperature, absolute temperature...) and, at the 
same time, we defined the ECMWF's analyses limitations: 
we concluded that in the first 10-20 m, the ECMWF analy- 
ses show a discrepancy with respect to measurements of the 
order of 2-3 m.s~^ for the wind speed and of 4-5 K for the 
temperature. 

The Meso-Nh model has been proven to be re- 
liable in re producing 3D maps of optical turbulence 
(jMasciadri et al. 1999a. b, 2001) and it has been statistically 
validated dMasciadri and Jabouilld I2OO1I : iMasciadri et al.l 
I2OO4I : iMasciadri and Egnerll2006h "alDOve mid-latitude astro- 
nomical site£]. Preliminary tests concerning the optimiza- 

^ See details in the next section 



tion of the model configuration and sensitivity to the hor- 
izontal and the verti cal resolution have a lready been con- 
ducted by our team ([Lascaux et al ][2QQ3) for the internal 
Antarctic plateau. In this paper we intend to quantify the 
performances of the model above this peculiar environment. 
More precisely, our goals are: 

- to compare the performances of the mesoscale Meso-Nh 
model and the ECMWF GCM in reconstructing wind speed 
and absolute temperature (main meteorological parameters 
from which the optical turbulence depends on) with respect 
to the measurements. This analysis will quantify the per- 
formances of the Meso-Nh model with respect to the GCM 
from the ECMWF. 

- to perform simulations of the optical turbulence above 
Dome C (75°06'04"S, 123°20'48"E) employing different 
model configurations and compare the typical simulated 
thickness of the surface layers well as the seeing in the free 
atmosphere with the one measured by Trinquet et al. (2008') 
(hereafter TR2008). In this way we aim to establish which 
configuration is necessary to reconstruct correctly the Cj^. 
In summary we aim to validate the Meso-Nh model on the 
Antarctic site. 

The two issues: (1) the surface layer thickness hsi and (2) 
the typical seeing in the free atmosphere are certainly the 
two main features that might get this place on the Earth 
extremely appealing for astronomers and it might be ex- 
tremely useful to have an independent confirmation from 
models of the typical values measured on the site. This study 
is focused on the winter season. In Section 2 we present the 
Meso-Nh model and the different configurations that were 
used to perform numerical weather simulations above the 
internal antarctic plateau. Section 3 is devoted to a sta- 
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tistical comparison of standard meteorological parameters 
(wind speed and temperature) deduced from Meso-Nh sim- 
ulations, ECMWF analyses and radiosoundings. In Section 
4 we present the results of the computation with Meso-Nh 
of the surface layer thickness for 15 nights in winter time 
and a comparison with the observed surface layer thickness 
from TR2008. Finally conclusions are drawn in Section 5. 



2 MODEL NUMERICAL SET-UP 

Meso-Nh (jLafore et alJ Il998l l is the non-hydrostatic 
mesoscale research model developed jointly by Meteo- France 
and Laboratoire d'Aerologie. 

It can simulate the temporal evolution of the three- 
dimensional atmospheric flow over any part the globe. The 
prognostic variables forecasted by this model are the three 
cartesian components of the wind v, w, the dry potential 
temperature G, the pressure P, the turbulent kinetic energy 
TKE. 

The system of equation is based upon an anelastic 
formulation allowing for an effective filter ing o f acous- 
tic waves. A Gal-Chen and SommervillJ (|l975l ) coordi- 
nate on the vertical and a C-grid in the formulation of 
Arakawa and Messingei (1976) for the spatial digitaliza- 
tion is used. The temporal scheme is an explici t three- 
time- level leap-frog scheme with a time filter (|Asselinl 
Il972l ). The turbulent scheme is a one-dimensional 1.5 clo- 
sure scheme (Cuxart et al. 2000) with the Bougeault and 
Lacarrer4 (|1989I ) mixing length. The surface exchanges are 
computed in an externalized surface scheme (SURFEX) 
i ncluding different physical packages, among which ISBA 
(|No ilhan and Planton 1989) for vegetation. 

Masciadri et al. (1999a,b) implemented the optical tur- 
bulence package to be able to forecast also the optical tur- 
bulence {C% 3D maps) and all the astroclimatic parameters 
deduced from the C%. We will refer to the 'Astro-Meso-Nh 
code' to indicate this package. To compare simulations with 
measurements the integrated astroclimatic parameters are 
calculated integrating the C% with respect to the zenith 
in the Astro-Meso-Nh code. The parameterization of the 
optical turbulence and the reliability of the Astro-Meso- 
Nh model have been proved in successive studies in which 
simulations have been compa red to measurement s pro- 
vided by different instrume nts ([Masciadri et al ] l200lL l2004l : 
iMas ciadri and Egner> l2006h . This has been achieved thank 
to a dedicated calibration procedur e that has been proposed 
and validated by the same authors (IMasciadri and Jabouillj 
[2001 ). 

The atmospheric Meso-Nh model is conceived for re- 
search development and for this reason is in constant evo- 
lution. One of the major advantages of Meso-Nh that was 
not available at the time of the Masciadri's studies is that it 
allows now for the use o f the interactive grid- nesting tech- 
nique (|Stein et al.l |2000[ ) . This technique consists in using 
different imbricated domains with increasing horizontal res- 
olutions with mesh-sizes that can reach 10 meters. 

We use in this study the Astro-Meso-Nh package, imple- 
mented in the most recent version of the atmospheric Meso- 
Nh model. To facilitate the put in the context of this work, 
the differences that have been implemented in the model 



configuration with respect to the previous Masciadri's stud- 
ies are listed here: 

• A higher vertical resolution near the ground has been 
selected. We still work with a logarithmic stretching near 
the ground up to 3.5 km but we start with a first grid point 
of 2 m (instead of 50 m) with 12 points in the first hundred 
meters. This configuration has been allowed thanks to the 
extremely smooth orography of this region of the Earth. It 
is obviously preferable because it permits to better quan- 
tify the turbulence contribution that typically develops in 
the thin vertical slabs in the first hundred of meters above 
the internal Antarctic Plateau. Above 3.5 km the vertical 
resolution is constant and equal to Aif=600 m as well as 
in Masciadri's previous work. The maximum altitude is 22 
kilometers. 

• The grid-nesting (see Table [1]) is implemented with 3 
imbricated domains allowing a maximum horizontal reso- 
lution of 1 km in a region around the Concordia Station 
(80 km X 80 km). 

• The simulations are forced at synoptic times (every 6 
hours) by analyses from the ECMWF. This permits to per- 
form a real forecast of the optical turbulence. To avoid mis- 
understandings, we highlight indeed that, as it has been 
extensively explained in previous studies (Masciadri et al. 
2004, Masciadri & Egner 2006), the Meso-Nh model has 
been used so far for simulations of the optical turbulence 
in a configuration permitting a quantification of the mean 
optical turbulence during a night and not a forecast of the 
optical turbulence. We perform therefore a step ahead with 
respect to results obtained so far with the Astro-Meso-Nh 
code. 

In spite of the fact that the orographic morphology is 
almost flat above Antarctica, it is known that even a weak 
slope can be an important factor to induce a change in the 
wind speed at the surface in these regions. The physics of 
the optical turbulence strongly depend on a delicate balance 
between the wind speed and temperature gradients. In or- 
der to study the sensitivity of the model to the horizontal 
resolution and to identify which configuration provides more 
reliable estimates we performed two sets of simulations with 
different model configurations. 

In the first configuration (that we will call monomodel) 
we used an horizontal resolution AX=100 km covering the 
whole Antarctic continent (Figure [T^,b and Table [1]) . We 
selected this configuration because it permits us to discuss, 
where it is possibl e, our r esults with respect to those ob- 
tained by Swain & Galled (2006) (hereafter SG2006) with 
the regional atmospheric model MAR above the Antarc- 
tic Plateau. In that case, indeed, the authors used this ex- 
tremely low horizontal resolution that has the advantage to 
be cheap from a computational point of view. This model 
configuration permits fast simulations but it is certainly nec- 
essary to verify that it is high enough to correctly resolve 
the most important features of the optical turbulence near 
the ground and in the high part of the atmosphere. 

In the second configuration (Table [1]) we used the grid- 
nesting technique, more expensive from a computational 
point of view but potentially more accurate in the recon- 
struction of the spatial distribution of the optical turbu- 
lence. The grid-nested simulations involved three domains. 
The largest domain covers all the Antarctic Plateau with 
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Table 1. Meso-Nh model configuration. In the second column the 
horizontal resolution AX, in the third column the number of grid 
points and in the fourth column the horizontal surface covered by 
the model domain. 





AX (km) 


Grid Points 


Surface (km) 


Monomodel 


100 


60x60 


6000x6000 




25 


120x120 


3000x3000 


Grid-Nesting 


5 


80x80 


400x400 




1 


80x80 


80x80 



120x120 points and it has a 25 km mesh-size (Figure [T]:). 
The second domain has a 5 km mesh-size, 80x80 points and 
is centered above the Dome C (Figure [TJi). The innermost 
domain has a 1 km mesh-size, 80x80 points and is centered 
above the Concordia Station area near the Dome C (Fig- 
ure [T^). Due to the fact that the typical topography of the 
internal Antarctic Plateau is much smoother than what is 
typically observed at mid-latitude sites, it has been decided 
to use a maximum resolution of 1 km instead of 500 m as it 
has been done in all the previously Masciadri's cited studies. 

The use of high-resolution has one first major impact: 
as it can be seen in Fig. [T] the Dome C area is more fairly 
reproduced in the grid-nested simulation than in the low 
horizontal resolution simulation (Figure [T]3,e). The altitude 
above mean sea level of the Concordia Station with high 
resolution is around 3230 m, whereas it is around 3200 m 
with the low resolution grid. 



3 MESO-NH SIMULATIONS: ABSOLUTE 
TEMPERATURE AND WIND SPEED 

The purpose of this section is to verify the performances of 
the mesoscale Meso-Nh model above the internal Antarctic 
Plateau and to verify if such a mesoscale model can provide 
a better estimate of the atmospheric flow than the GCM 
from the ECMWF. 

An important number of winter nights (47) were simu- 
lated with the Meso-Nh mesoscale model. We analyze here 
the key meteorologic parameters from which the optical tur- 
bulence depends on: the temperature and the wind speed. 
Both configurations (low horizontal resolution monomodel, 
and high horizontal resolution grid-nesting) are tested and 
evaluated. The model is initialized with ECMWF analyses 
extracted in the nearest grid point with respect to Dome 

All the simulations started at 00 UTC and were inte- 
grated for 12 hours. Simulations outputs at 12 UTC are 
compared with measurements we retrieved from the site 
\http ://www. dimantartide. it) as well as with the analyses 
from the GCM of the ECMWF. Every 6 hours we forced 
the simulations with the ECMWF analyses in order to avoid 
that the model diverges and/or correct the atmospheric flow 
as a function of the predictions at larger spatial scales. 

In this section a statistical study of the wind and tem- 
perature profiles at Concordia Station, Dome C, is per- 
formed. The 47 nights have been selected in June, July and 

^ Horizontal resolution of the ECMWF analyses: ~0.5° 



Table 2. Mean values on 47 winter days at Concordia Station 
near Dome C of wind speed and temperature at the surface 
level, for radiosoundings and ECMWF analyses. Into brackets, 
the corresponding statistical error (ct/a/TV). 





Radiosondes 


ECMWF 


Wind speed (m.s ^) 
Temperature (K) 


4.02 (± 0.37) 
212.90 (± 1.11) 


6.51 (± 0.37) 
216.64 (± 0.85) 



Table 3. Mean values on 47 days at Concordia Station near 
Dome C of wind speed and the temperature at the surface 
level, for the Meso-Nh simulations: 1-MOD is for the low hor- 
izontal simulation with AX=100 km and Grid-N is for the 
high horizontal grid-nested simulation with AX=1 km for the 
innermost model centered above the Dome C area. Into brack- 
ets, the corresponding statistical error (a/VN). 





1-MOD 


Grid-N 


Wind speed (m.s ^) 


4.23 (±0.26) 


3.98 (± 0.28) 


Temperature (K) 


214.92 (± 0.68) 


214.50 (± 0.72) 



August 2005 and July 2006. For all the 47 nights selected, 
we respected the following criterion: 

• A radiosounding is available at the end of the simulation 
(at 12 UTC of the selected night) to perform comparisons 
between Meson-Nh outputs, ECMWF analyses and obser- 
vations. 

• We selected nights in which, the corresponding ra- 
diosoundings launched at time t=to+12 h (with to the initial 
time of the simulation) cover the longest path along the z- 
axis (perpendicular to the ground) before to explode. It was 
impossible to collect in winter time 47 nights in which all 
the balloons reached 20 km. The mean altitude reached by 
the selected balloons was about 10 km above ground level. 

3.1 Model validation: vertical profiles of 
temperature and wind speed 

Figures [2] and [3] show the mean vertical profiles of temper- 
ature and wind speed, respectively, computed for 47 winter 
nights from the two model configurations (low and horizon- 
tal resolution), the ECMWF analyses and the radiosound- 
ings. The location of the profiles is Concordia Station, in 
the Dome C area. All profiles have been interpolated on a 
regular 5 m vertical grid, in order to ease the comparison. 

3.1.1 Temperature 

The mean temperature profiles are very similar over the en- 
tire free atmosphere (Figure [2^). In the first kilometer the 
temperature gradients reconstructed by the Meso-Nh sim- 
ulations (with high and low resolution) and the ECMWF 
analyses are not as pronounced as the one obtained with 
the radiosoundings. This means that the mesoscale model 
as well as the General Circulation Model (ECMWF) recon- 
struct a shghtly less stable atmosphere in this region even 
if the mesoscale model better approaches the observations 
trend. However, as it will be shown in the next section, the 
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Figure 2. Mean temperature profiles at the Concordia Station in the Dome C area over 47 winter nights (bold lines) and the corresponding 
standard deviation (dashed lines). Radiosoundings: black lines; ECMWF analyses: blue lines; Meso-Nh with low horizontal resolution 
(monomodel): green lines; Meso-Nh with high resolution (grid-nesting): red lines. Profiles are displayed up to (a) 20 km, (b) 1 km and 
(c) 250 m above ground level. Units in Kelvin (K). 
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Figure 3. Mean wind speed profiles at the Concordia Station in the Dome C area over 47 winter nights (bold lines) and the corresponding 
standard deviation (dashed lines). Radiosoundings: black lines; ECMWF analyses: blue lines; Meso-Nh with low horizontal resolution 
(monomodel): green lines; Meso-Nh with high-resolution (grid-nesting): red lines. Profiles are displayed up to (a) 20 km, (b) 1 km and 
(c) 250 m above ground level. Units in m.s"-*^. 



mesoscale model (Meso-Nh) provides a more accurate esti- 
mate of the surface temperature than what the ECMWF 
model can do. 

3.1.2 Wind speed 

Figure [3] shows the mean wind speed during the 47 days, 
with the corresponding standard deviation. From the ground 
up to 10 km, analyses and radiosoundings mean wind speed 
are well correlated. Above 10 km the wind speed recon- 
structed by the ECMWF analyses is slightly larger than the 
one reproduced by Meso-Nh (monomodel and grid-nesting 
- left of figure El). It is hard to say whether the ECMWF 
analysis or the Meso-Nh simulation is the best since no 
mean value from the observations is available at this alti- 
tude. Between 1 km and 10 km there are no major differ- 
ences between the mesoscale model and the ECMWF. Below 
1 km it is well visible that Meso-Nh better reconstructs the 
strong wind shear than the ECMWF analyses up to achieve 
a more correlated wind speed value near the surface. At 
150 m this difference is maximized: the wind speed provided 
by the ECMWF analyses is a bit too weak (12 m.s~^ instead 
of 14 m.s~^ in the observations). At the same altitude the 
Meso-Nh simulations give better results, with a mean wind 
profile perfectly correlated to the one measured by the ra- 
diosoundings. The improvement is even better in the case of 



the high-resolution model. The difference between low hor- 
izontal and high horizontal simulations is more important 
above 12 km, with an increase in intensity of the wind more 
important in the high-resolution simulation. These results 
match pe rfectly with a dedicated analysis that our group 
did (HageUn et al.ll2008l l on a comparison between the wind 
speed provided by ECMWF analyses and radiosoundings 
near the surface. 



3.2 Model validation: the surface 

The mean values of the surface wind speed and absolute tem- 
perature at Dome C were computed for the 47 nights. The 
results are reported in the Tables [2] and [3] (mean values for 
ECMWF analyses, radiosoundings, and Meso-Nh low and 
high horizontal res olution simul ations, respectively). In an- 
other paper ( Hagelin et al.ll2008l l our group showed that the 
radiosoundings in the first grid-point perfectly match with 
measurements provided by the automatic weather stations 
(AWS). 

3. 2. 1 Temperature 

One can see that the ECMWF an alyses, as already rep orted 
in a previous paper of our team ([Hagelin et al.|[2008f ) on a 
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different data-set sample, are too warm at the surface (with 
a difference of almost 4 K, table [2]) in winter with respect to 
the observations. 

The mean surface temperature simulated by Meso-Nh 
after 12 hours is closer to the observations than the ECMWF 
analyses. The difference between ECMWF analyses and 
observed mean temperature is ATecmt(;/,obs=3.74 K. The 
mean surface temperature in the low-resolution simulation is 
ATmn^-Zot(;,obs=2.02 K higher than in the observations. It is 
only ATmn^-^ip^,o6s=l-60 K higher for the grid nested sim- 
ulation. This means that the mesoscale model reconstructs 
a surface temperature that is typically a factor ^^2-2.5 (that 
is /Sr£ec7nwf,ohs/^^7nnh-high,ohs) morc accuratc than the 
ECMWF analyses. More over it can be seen that the sur- 
face temperature is better retrieved with the use of the high 
horizontal resolution [/\X—1 km in the innermost domain) 
than the low horizontal resolution (AX=100 km). 



a 



High 
Low 



50 100 

Observed Surface Layer Thickness 



150 



3.2.2 Wind speed 

Tables [2] and [3] show the mean values of the wind speed 
at the first interpolated point of the profiles. The mean 
wind speed in the ECMWF analyses (6.51 m.s~^) is higher 
than the observed wind speed (4.02 m.s~^), thus a dif- 
ference of [/^Y ecmw f ,ohs— 2.49 m.s~^jf|. Both mesoscale 
low and high horizontal simulations reproduce more accu- 
rately the surface wind speed than the ECMWF analy- 
ses. With a mesh-size of 100 km in Meso-Nh, the differ- 
ence between the simulated and measured mean wind speed 
is of /SN mnh-iow ,ohs— 0.21 m.s~^. The grid nested simula- 
tions [/\X—1 km in the innermost domain) give even bet- 
ter results with a difference of /SN mnh-high,ohs—^-^^ m.s""*^ 
only. This means that the mesoscale model reconstructs a 
surface wind speed that is typically a factor ~ 60 (that 
is ecmwf,ohs/ ^Ymnh-high,ohs) morc accuratc than the 
ECMWF analyses. 

The median value obtained with the mesoscale model 
MAR (SG2006) differs from the measurements, done on 
the same statistical sample, /SN MAR,ohs— 0.9 m.s""*^. We 
conclude therefore that our simulations performed with 
Meso-Nh are a factor ~ 4 {/^Y MAR,ohs/ mnh-iow,ohs) 
more accurate than those performed with MAR by SG2006 
if we use the same horizontal resolution. They are a factor 
~ 22 {AVMAR,obs/^'Vmnh-high,obs) morc accuratc if we 
use the high horizontal resolution with Meso-Nh. We also 
note that in that paper is reported just the standard 
deviation (a) and not the statistical error {cf/^/N) where 
N is the number of independent estimates. The former is 
appropriate to describe the dispersion of the data-set and 
the latter is useful to describe the precision of the estimate 
of the mean values and it is therefore more appropriate 
to quantify the performances of the simulations with 
respect to observations. This means that each standard 
deviation in Table 1 (SG2006) should be multiplied by 



Difference from the iHagel m et al.l (|2Q08l ) paper values are just 
due to the fact that in that paper, all the nights of the three 
months (June, July and August) were used while in this paper 
we simulated just 47 nights selected in two different years (2005 
and 2006). The statistical sample is therefore not the same. 



Figure 4. Correlation plot between measured and simulated sur- 
face layer thicknesses (black: monomodel configuration; red: grid- 
nested configuration). For the simulated values only the mean 
values between 12 UTC and 16 UTC are considered. For each 
configuration of the simulation (high and low horizontal resolu- 
tion) the error bars are reported for one point only (and are equal 
to a). Units are in meter (m). 

1 / a/IV, where N=90 is the number of nights in three months. 



4 MESO-NH SIMULATIONS: OPTICAL 
TURBULENCE 

We have seen in the previous section that the Meso-Nh 
model allows for a better forecasting of meteorological pa- 
rameters such as wind speed and temperature than the anal- 
yses from the GCM of the ECMWF. These thermodynamic 
parameters are very important for the forecasting of optical 
turbulence because the computation of the astroclimatolog- 
ical parameters (seeing e, isoplanatic angle ^o, wavefront co- 
herence time To) depend directly from them. As it has been 
highlighted in the introduction, the model has been run with 
the Astro-Meso-Nh package that allows the prediction of 3D 
maps. The most important features that characterize 
the optical turbulence above the internal Antarctic Plateau 
are: 

• the typical surface layer thickness 

• the median seeing in the free atmosphere 

• the median seeing in the whole atmosphere. 

These three elements permit us to do a complete analysis of 
the optical turbulence developed above Dome C. Our tests 
are made on all the 15 nights in winter time for which mea- 
surements of the optical turbulence surface thickness, and 
partial seeing in different vertical slabs (free atmosphere see- 
ing and total seeing) are available (see TR2008). Due to the 
fact that in the Trinquet et al. (2008) paper is available just 
the median profile for the four seasons and not the Clj 
profile for each individual nights, we selected our sample 
taking all and only the nights belonging to the winter as 
defined by TR2008 i.e. between June 21 to September 21. 
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Table 4. Surface layer thickness hgi for 15 winter nights 
(TR2008). Units in meter (m). The criterion used is the one from 
Eq[T] The mean value is also reported with the associated statisti- 
cal error ct/^/N where N is the number of independent estimates. 



Date 


Observed surface 


Date 


Observed surface 




layer thickness 




layer thickness 


04/07/05 


30 


12/08/05 


22 


07/07/05 


21 


29/08/05 


47 


11/07/05 


98 


02/09/05 


41 


18/07/05 


26 


05/09/05 


20 


21/07/05 


47 


07/09/05 


39 


25/07/05 


22 


16/09/05 


24 


01/08/05 


40 


21/09/05 


22 


08/08/05 


30 








Mean 


35.3 






(7 


19.9 








5.1 





This is the most interesting period for stellar astronomical 
applications. 



4.1 Surface layer thickness (hsi) 

In order to verify how well the simulated surface thickness 
(hsi) matches with the measured one we are forced to com- 
pute the typical height of the surface layer using the same 
criterion as in TR2008. The authors defined the thickness 
hsi of the surface layer as the part containing 90% of the 
first kilometer boundary layer optical turbulence: 



< 0.90 



(1) 



where C"^ is the refractive index structure parameter. 

The observed hsi for 15 winter nights from TR2008 
are reported on table IH The simulated hsi for the same 
15 nights, calculated for the two configurations (monomodel 
i.e. low horizontal resolution and grid-nesting i.e. high hori- 
zontal resolution), are reported in Table [5] 

Two different time intervals were chosen for the com- 
putations of the mean values (11 UTC to 17 UTC and 12 
UTC to 16 UTC), both temporally centered around 14 UTC 
i.e the time at which the balloons have been launched. We 
reported on Figure |4] a correlation plot where the measured 
and the simulated (mean values between 12 UTC and 16 
UTC only) hsi are compared. 

It is worth to highlight how the comparison between 
measurements and simulations is done. At present time, it 
is meaningless to predict the C'^ at a time t=to. The opti- 
cal turbulence is indeed a parameter fiuctuating in space and 
time at much smaller scales than what a classical meteoro- 
logical parameter does and at much smaller scales than the 
model mesh sizes. The optical turbulence is parameterized in 
the model and not resolved explicitly and it has to be quan- 
tified statistically. For this reason it is not realistic at present 
time, to forecast the Cjj profiles with a better precision in 



Table 5. Mean surface layer thickness hsi for the same 15 winter 
nights than in tablelH deduced from Meso-Nh computations using 
the criterion in Eq. ^ Two different times intervals were chosen: 
between 11 UTC and 18 UTC, and between 12 UTC and 16 
UTC. Units in meter (m). The mean value is also reported with 
the associated statistical error cr/yfN. 



Date 


Surface layer thickness 


Surface layer 


thickness 




Meso-Nh grid-nesting 


Meso-Nh monomodel 




11-17 UTC 


12-16 UTC 


11-17 UTC 


12-16 UTC 


04/07/05 


24.5 


25.3 


43.9 


44.9 


07/07/05 


37 


36.9 


33 


32.9 


11/07/05 


65.1 


63.9 


91.7 


91.7 


18/07/05 


38.9 


38.2 


72.4 


72.3 


21/07/05 


46.4 


46.2 


56.9 


56.8 


25/07/05 


83.4 


82.1 


31.6 


32 


01/08/05 


27.9 


27.9 


41.1 


41 


08/08/05 


33.3 


29.7 


145.7 


153.3 


12/08/05 


23.1 


18.2 


44.3 


44.8 


29/08/05 


92 


84.6 


77.3 


78 


02/09/05 


52.8 


54.5 


66.5 


67.1 


05/09/05 


120 


122 


37.2 


38.8 


07/09/05 


58.1 


56.9 


89 


85.9 


16/09/05 


27.6 


29.4 


40 


41.7 


21/09/05 


17 


17.3 


105 


107.5 


Mean 


49.8 


48.9 


65 


65.9 


(7 


29.4 


29.3 


32.7 


33.6 


a/VN 


7.6 


7.6 


8.2 


8.7 


time than 


a At of a 


few hourfl Our objective 


is to find 



therefore a correlation between the measurements obtained 
during one balloon launch (C^ profile) and the mean of 
the temporal evolution of the simulations (C^ profiles) ex- 
tended on a At of a few hours. Looking at Table [5] we observe 
that the values of the simulated surface layer thickness are 
weakly dependent on the selection of the temporal interval of 
integration (At: 11-17 UTC and 12-16 UTC). We chose the 
second one because the associated simulations give a slightly 
better correlation with the measurements (Table 2]). Using 
the criterion expressed by Eq. [1] the grid-nested simulations 
give a mean surface thickness hsi,mnh-high— 48.9±7.6 m 
(where (t/^/N—7.Q m) and the monomodel a mean thick- 
ness hsi,mnh-iow= 65.9=b8.7 m (Table [5]). The low resolu- 
tion configuration leads to a higher mean forecasted thick- 
ness with respect to the observed one (/isi,obs =35. 3±5.1 m) 
while the grid-nested technique is closer to the observations 
{Ahobs,mnh-high =13.6 m) and within the typical a. If we 
take into account the statistical error m) we 

conclude that the mesoscale model provides, for this statisti- 
cal sample of 15 nights, a typical surface layer thickness just 
~ 6 m higher than the observed one. The dispersion a of the 
surface layer thickness for the simulated nights (cr=29.3 m) 
is just slightly larger than the observed one (cr=19.9 m). This 
indicates that the intrinsic dispersion of the hsi is well recon- 
structed by the model. Fig. 2] shows the correlation between 
simulated and observed hsi with the corresponding a values. 



^ We refer the reader to iMasciadri and Egnerl (|2p06')-Session 2 
for an extended explanation of this concept. The perspective of 
this typology of studies is to attain smaller and smaller At. 
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Figure 5. Temporal evolution of the profile at Concordia Station, Dome C obtained with Meso-Nh, between 00 UTC and 18 UTC for 
the three nights: 4/7/2005, 18/7/2005 and 12/8/2005 (from the top to the bottom). On the left-side are shown the simulations obtained 
with the low horizontal resolution, on the right-side the simulations obtained with the high horizontal resolution. The mean value of the 
surface layer computed with the criterion ^ between 12 UTC and 16 UTC is reported for each night. The thin black line represents the 
evolution of the surface layer height during the night. The first couple of hours can be considered as spurious values because of the model 
adaptation to the ground. 



We observe moreover, that, in spite of the more expensive 
simulations in terms of computing resources (time and mem- 
ory), the high horizontal resolution grid- nested configuration 
seems to be necessary to better reconstruct the concentra- 
tion of the turbulence in a thin layer near the surface. More 
precisely a horizontal resolution of 100 km provides a bias 
in the typical hsi of the order of ^ 30 m. This result is rele- 
vant with respect to the study done by SG2006 who used a 
resolution of 100 km and found a mean 39 m on a sample of 
90 nights. Unfortunately, the authors define the hsi as the 
elevation (starting from the lowest model layer) at which 
the turbulent kinetic energy contains 1% of the turbulent 
kinetic energy of the lowest model layer. This definition is 



completely different from that used by TR2008 (and by us) 
i.e. the elevation in which is included 90% of the optical 
turbulence developed in the first kilometer. It is therefore a 
no sense whatever comparison of measured hsi with the hsi 
by SG2006. To have an idea of the impact of the selected 
criteria on the mean hsi estimate we repeated the statistical 
calculation using the same criterion employed by SG2006. 
The mean hsi obtained with the high resolution is 32.3 m 
that is ^ 17 m lower than the 48.9 m estimated with the 
criterion used by TR2008. We therefore deduce that the typ- 
ical thickness estimated by SG2006 using the same criterion 
employed by TR2008 should be certainly much larger than 
39 m. This confirm the fact that mesoscale models provides 
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Figure 6. Left: Median profile measured (black line) with microthermal sensors mounted on balloons (from Trinquet et al. 2008) 
and simulated with the Meso-Nh mesoscale model with the low-horizontal resolution (blue line) and the high-horizontal resolution (red 
line) Right: zoom in the first 100 m. Units in m.~^/^. 



an overestimation of the hsi if they are run with a AX= 100 
km. 

What about the morphology of the surface layer ? Looking 
at Table 2] it is well visible that the dispersion (a) of the 
measurements is relatively large and sometimes the surface 
layer thickness can be of the order of many tens of meters. 
Is the mesoscale model Meso-Nh able to reconstruct such 
a large variability in the morphology of the surface layer ? 
Table [5] shows that basically for all the nights, the model 
reconstructs within the a value the corresponding observed 
value. In just one case (5/9/2005) the simulated thickness 
is much larger {r^ 120 m instead of the observed 20 m). We 
highlight that a larger surface thickness does not mean nec- 
essarily a thicker and developed turbulence near the ground 
but simply that in the (0, h^;) range there is 90% of the 
turbulence developed in the first kilometer. 

To highlight the potentiality of Meso-Nh in discriminat- 
ing between different turbulent nights, and to better visu- 
alize the impact of the model configuration (low and high 
resolution) on the forecasted optical turbulence, the tem- 
poral evolution of the profiles near the surface and the 
corresponding surface layer heights are displayed for three 
nights (Fig. O. The first hours of the simulation should not 
be taken into account because they correspond to the spin- 
up of the model, and do not describe realistic C% profiles. 
Fig. O left-side shows the simulations with low resolution, 
Fig.[5l right-side shows the simulations with high resolution. 
The three selected nights are, from the top to the bottom: 
4/7/2005, 18/7/2005 and 12/8/2005. In all the three cases 
the developed turbulence layer is thinner for the high hori- 
zontal resolution case than for the low resolution one. The 
high resolution provides also a better correlation to obser- 
vations. 

• During the night 4/7/2005 (Fig. [5^,b for monomodel 
and grid- nested configurations, respectively) the predicted 
surface layer remains constant in time, with a mean 
forecasted thickness well correlated to the observed one 



sl ,mnh — high 



^25.3 m, h. 



sl ,mnh — low 



^44.9 m and h 



sl,obs- 



30 m) . The high resolution simulation matches the observa- 
tions within a/VN. 

• During the night 18/7/2005 (Fig. [5t,d for monomodel 
and grid-nested configurations, respectively) the predicted 
surface layer remains constant in time. The simulation at low 
horizontal resolution {hsi,mnh-iow —72.3 m) overestimates 
the observation {hsi^ohs— 26 m) while the simulation with 
high resolution {hsi,mnh-high — 38.2m) is better correlated. 

• During the night 12/8/2005 (Fig. [5^,f for monomodel 
and grid-nested configurations, respectively) it is visible that 
the morphology of the reconstructed surface layer is thin as 
well as in the previous cases. For this night, the high horizon- 
tal resolution permits to put in evidence a better sensitivity 
of the model to the temporal variability of the that 
has been observed also in other nights. Indeed we observe, 
just above the surface layer in the last part of the simulation, 
short-bumps of optical turbulence i.e. fluctuations of the Cfj 
profile forecasted by the model near the surface. This 
variability is the signature of an evident temporal evolu- 
tion of the turbulent energy distribution even in conditions 
of a strongly stratified atmosphere. Also in this last case, 
even if both simulated heights {hsi^mnh-high —18.2 m and 
hsi,mnh-iow —44.8 m) are well below 100 m, the high reso- 
lution provides a clearly better correlation (within (t/\/N) 
to the observations {hsi,obs= 22 m). 



4.2 Optical turbulence vertical distribution: 
seeing in the free atmosphere and in the 
whole atmosphere 

Figure [6] shows the median of the C% profile measured 
by the microthermal sensors mounted on the balloons (15) 
launched at Dome C during the winter 2005 and simulated 
by the Meso-Nh model with the low and high horizontal 
resolutions. We observe that the shape of the C% is well 
reconstructed all along the 13 km by the model. Also the 
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model can reconstruct the above 13 km. In this region 
a comparison with measurements is not possible because the 
balloons usua ll y exp lode at these heights (see discussion in 
iHagelin et alJ ( 2008| )). In the first kilometer, the simulated 



CIj profile decreases (from the ground to higher altitudes) 
in a less sharp way than what monitored by observations 
(at 30 m from the ground). This is not surprising and it 
derives from the fact that simulations are slightly less ther- 
mally stable than the observations near the surface. How- 
ever, in the zoom of the first 100 m (Fig. (6] right-side) it is 
well visible that the shape of the median simulated pro- 
file done with the high horizontal resolution configuration is 
much better correlated with the shape of the observed me- 
dian profile than the one obtained with a low horizontal 
resolution. 

The seeing in the free atmosphere and in the whole atmo- 
sphere for A=0.5xl0~^m are: 



Cl{h)-dh 



3/5 



£tot — 5.41 • A 



\Jsra 



Cl{h)'dh 



3/5 



(2) 
(3) 



with htop ^ 13 km from the sea level i.e. where the balloons 
explode and we have no more their signal. 

Table [6] reports, for the 15 nights, the seeing in the 
free atmosphere {eFA=S[hsiMop]) whole atmo- 

sphere {£TOT=S[srn,htop]) Calculated for the simulations and 
the observations. Values of the observed stot are taken 
from TR2008. We considered \isi= 33m instead of 35.3 m 
(Table |4]) for the observations because, in the Trinquet et 
al. (2008) paper, the authors provide the seeing integrated 
above h^^, where h^^ is calculated on all the nights belong- 
ing to the [February- November] range. We considered for 
the simulations h^^ as retrieved from Table [5] 

Again we observed that results are weakly dependent 
on the temporal range on which the means values are cal- 
culated and for this reason we report just the 12-16 UTC 
case. Fig. [71 shows the correlation between the observed and 
simulated values for the seeing in the free atmosphere and 
in the whole atmosphere. The median of the observed see- 
ing in the free atmosphere for the 15 nights is epA^ohs — 
0.3d=0.2 arcsec; the median seeing in the free atmosphere 
simulated by Meso-Nh with the high horizontal resolution is 
h-hiah^ 0.35it0.24 arcscc and with the low horizon- 
tal resolution is eFA,mnh-high= 0.42zb0.28 arcsec. Both me- 
dian simulated values (with low and high resolution) match 
the median value obtained with observation within the sta- 
tistical error even if the high resolution is much better cor- 
related (relative error of 16%, Asobs,sim=0.05"). If we look 
at the total seeing developed on the whole atmosphere it 
is well visible (Table [6] and Fig]?]) that the model overes- 
timates the measurements with both resolutions. We have 
a simulated median STOT^mnh-iow— 3.58=b0.42 arcsec and 
£TOT,mnh-high= 2.29±0.38 arcscc versus an observed me- 
dian £TOT,ohs— 1.6i:0.2 arcsec. Even if we take into account 
the more accurate estimates (high resolution) we obtain a 
dispersion simulations/observations As ^ 0.7 arcsec. The 
excess of optical turbulence reconstructed by the Meso-Nh 
model is clearly concentrated in the surface layer. We can 
not exclude an underestimate from measurements but there 
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Observed Total Seeing 
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Observed Free Atmosphere Seeing 

Figure 7. Correlation plot between measured and simulated total 
(top) and free atmosphere (bottom) seeing (black: monomodel 
configuration; red: grid-nested configuration). For the simulated 
values only the mean values between 12 UTC and 16 UTC are 
considered. For each configuration of the simulation (high and low 
horizontal resolution) the error bars are reported for one point 
only (and are equal to cr). Units are in arcsec. 



are, at present time, no major elements that lead to this 
assumption. We are working, on the contrary, on a paper 
to explain this discrepancy and overcome this limitation. 
Considering that we proved that the meteorological param- 
eters are well reconstructed by the Meso-Nh model near the 
surface (Section 3) and that the surface numerical scheme 
(Interaction Soil Biosphere Atmosphere ISBA) responsible 
of the control of the budget of the turbulent ground/air 
fluxes has been recently optimized for Antarctic applications 
(|Le Moigne et al.l (|2008l )) in the context of our project, we 
concentrated our attention on the dynamical and optical nu- 
merical turbulence schemes. 

In terms of comparison with the SG2006 study we note 
that the latter study indicates a typical underestimated to- 
tal seeing of 1.16 arcsec with respect to the observed one 
(1.6 arcsec). The discrepancy is smaller from a quantitative 
point of view (A£tot=0.45 arcsec) with respect to what 
we find and it is in the opposite direction. The questionable 
issue in the SG2006 study is that the turbulence kinetic 
energy provided by SG2006 in the flrst lev els of t h e MA R 
model is often of the order of 10"^ m^s-^ (jGalled (|2007f )). 
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Table 6. Total seeing eTOP=^[Sm,htop] seeing in the free 
atmosphere eFA=^[hsi,htop] calculated for the 15 nights and av- 
eraged in the temporal range 12-16 UTC. See the text for the 
definition of h^^ and htop- In the second column are reported the 
observed values, in the third and fourth columns the simulated 
values obtained with high and low horizontal resolution respec- 
tively. Units in arcsec. 





Obs. 


MESO-NH 


MESO-NH 






HIGH 


LOW 


Date 


sfa/^tot 


efaI^tot 


efaI ^tot 






(^sZ=48.9m) 




04/07/05 


0.3 / 1.6 


0.20 / 3.40 


0.45 / 4.61 


07/07/05 


0.2 / 1.5 


0.35 / 3.03 


0.31 / 2.64 


11/07/05 


1.4 / 1.7 


2.42 / 3.55 


3.27 / 4.32 


18/07/05 


0.3 / 2.0 


0.35 / 3.67 


1.57 / 4.68 


21/07/05 


0.7 / 1.1 


0.51 / 2.12 


0.52 / 2.17 


25/07/05 


0.3 / 1.0 


0.29 / 0.84 


0.35 / 3.58 


01/08/05 


0.8 / 1.6 


0.26 / 3.91 


0.28 / 4.17 


08/08/05 


0.5 / 2.3 


0.46 / 2.21 


0.46 / 1.12 


12/08/05 


0.2 / 1.5 


0.32 / 1.96 


0.42 / 5.16 


29/08/05 


2.5 / 3.6 


2.67 / 3.35 


1.97 / 4.52 


02/09/05 


0.9 / 1.9 


0.66 / 2.26 


0.73 / 2.82 


05/09/05 


0.3 / 1.0 


0.38 / 1.05 


0.37 / 2.81 


07/09/05 


1.4 / 2.8 


2.58 / 6.77 


3.33 / 7.02 


16/09/05 


0.2 / 1.5 


0.28 / 1.44 


0.27 / 2.19 


21/09/05 


0.2 / 1.7 


0.32 / 2.29 


0.27 / 0.87 


Median 


0.3 / 1.6 


0.35 / 2.29 


0.42 / 3.58 


(7 


0.7 / 0.7 


0.92 / 1.46 


1.07 / 1.64 




0.2 / 0.2 


0.24 / 0.38 


0.28 / 0.42 



This values is extremely low and it basically indicates no 
turbulent kinetic energy on the first level of the model and 
such a condition is contrary to what observed with measure- 
ments. This is consistent with the fact that the MAR model 
underestimates the seeing in the surface layer. 

We conclude that the Meso-Nh model, in the present 
configuration, reconstructs with good statistical reliability 
the hsi and the seeing in the free atmosphere while shows a 
tendency in overestimating the strength of the seeing in the 
surface layer. The interesting result of this paper is therefore 
the fact that the most important features for astronomical 
interest (the surface layer thickness and the typical seeing in 
the free atmosphere) observed with measurements are con- 
firmed with mesoscale atmospherical model. We note that 
the this is the first confirmation made by a mesoscale model 
of the typical seeing in the free atmosphere. Besides it is 
worth to highlight that these are the first C% simulations 
ever done above the internal Antarctic Plateau and extended 
all along the whole atmosphere. Figure [8] shows the tempo- 
ral evolution of the C% profile in the free atmosphere (more 
precisely in the (1,12) km vertical slab) related to three se- 
lected nights in the sample of the 15 simulated nights. In all 
of the three nights it is well visible that, even at such high 
altitudes, the model is active and the vertical distribution of 
the optical turbulence changes in time with a not negligible 
dynamic from a quantitative point of view. The Cfj values 
extend, indeed, on the logarithmic scale (-18,-16.5). In all 
the 3 cases it appears clearly that the high-horizontal res- 
olution provides a better temporal variability as expected. 



These results are therefore very promising in terms of pre- 
dictions of the Cj^ 3D maps on long time scales. 



5 CONCLUSIONS 

In this paper we study the performances of the Meso-Nh 
mesoscale meteorological model in reconstructing meteoro- 
logical parameters (wind speed and temperatures) as well 
as the optical turbulence above Concordia Station in the 
Dome C area, a site in the internal Antarctic Plateau. This 
is, at our knowledge, the first study concerning the optical 
turbulence reconstructed with an atmospherical mesoscale 
model above Antarctica on the whole atmosphere. This 
study is concentrated on the winter season i.e. the most 
interesting for stellar astronomical applications. The vali- 
dation of the model for the meteorological parameters has 
been done comparing measurements (radiosoundings) and 
simulations on a sample of 47 nights. The validation of the 
model for the optical turbulence has been done comparing 
simulations with measurements on a sample of 15 nights. 
Two different model configurations were tested: monomodel 
simulations using a low horizontal resolution (AX =100 km) 
and grid-nesting simulations with high horizontal resolution 
(AX=1 km for the innermost domain). The low resolu- 
tion model permitted us to discuss the results obtained 
previously in the literature. The observations used for 
the validation are, for the meteorological parameters, the 
analyses from the ECMWF Global Circulation Model and 
radiosoundings (47 nights) and, for the optical turbulence, 
the Cfj and seeing values (15 nights) measured in situ 
(Trinquet et al. 2008). 

The main conclusions of this study are: 

(1) We showed that near the surface, Meso-Nh retrieved 
better wind speed vertical gradient (wind shear) than the 
ECMWF analyses from a qualitative as well as quantitative 
point of view, thanks to the use of a highest vertical 
resolution. We expect therefore a better reconstruction 
of the katabatic winds typical of these regions by the 
Meso-Nh model than the GCM models. Also Meso-Nh 
better reconstructs the thermal stability near the sur- 
face than the GCMs. The analysis of the first vertical 
grid point permits us to conclude that the Meso-Nh 
model surface temperature is closest to the observations 
ATmnh-high,obs^l.QO K than the ECMWF General Cir- 
culation Model (ATecmt(;/,o6s=3.74 K) which is too warm. 
The improvement for the estimate of the wind speed is 
even more evident {/SN mnh-high,ohs— 0.04 m.s~"^ versus 
AVecmio/,obs= 2.49 m.s~^). 

(2) For what concerns the parameters concerning the 
optical turbulence, again the results are resolution de- 
pendent. The simulations with low resolution provides a 
too thick surface layer (almost double of the observed 
one) while those with high resolution provide a mean 
hsZ,mn^-^i^^=48.9±7.6 m versus an equivalent observed 
hsZ,o6s=35.3±5.1 m. Taking into account the statistical error 
we observe that the high horizontal mode provides a surface 
layer thickness that is statistically just 6 m higher than the 
observed one but within the dispersion a of the observations. 
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MONOMODEL GRID-NESTING 




Figure 8. Temporal evolution on 18 hours of the profiles in the vertical slab (1-12) km related to three nights chosen in the selected 
sample of 15 nights. The dynamic of the (color palette) has been tuned to put in evidence the variation of the values (in 

logarithmic scale) during the time in the range (-18,-16.5). 



(3) The integral of the above the hsi i.e. the seeing 
in the free atmosphere £FA,ohs= 0.3±0.20 arcsec is recon- 
structed with an excellent level of reliability (As=0.05 
arcsec) by the model used with the high resolution configu- 
ration eFA,mnh-high= 0.35±0.24 arcscc. The low resolution 
provides a worse estimate even if within the a of the 
observations. 

(4) The model still shows a tendency in overestimating the 
turbulence in the surface layer. For an observed STOT^ohs — 
1.6 arcsec we have a simulated £TOT,mnh— 2.29 arcsec with 
the model in high horizontal resolution mode. This is the 
subject of an on-going study conceived to answer to this 
open question. 



(5) The results concerning the computation of the mean 
thickness of the surface layer as well as the seeing in 
different vertical slabs are not very dependent of the time 
interval used to average it. This widely simplifies the 
analysis of simulations. 

(6) Estimates obtained with the grid-nested simulations 
are closer to the observations than those obtained with 
monomodel simulations. This study highlighted the neces- 
sity of the use of high horizontal resolution to reconstruct 
a good meteorological field as well as the parameters char- 
acterizing the optical turbulence in Antarctica, even if the 
orography is almost flat over the internal Antarctic Plateau. 
The employment of the low resolution (100 km) alone can 
hardly be used to identify the best site on the Antarctic 
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Plateau. However, it can be used to identify rapidly, on the 
whole Antarctic Plateau, the most interesting regions in 
which to focus, successively, simulations at high horizontal 
resolutions on smaller surfaces domains. With "the most 
interesting regions" we mean those with the lowest surface 
layer thickness for example. 

(7) The Meso-Nh model is able to reconstruct a mean 
C]^ profile well fitting the vertical optical turbulence 
distribution measured in the first 20 km from the ground. 
The model also shows a not negligible temporal variability 
in the whole 20 km from the ground in a very small dynamic 
range. The latter is to be considered a very interesting 
feature because it is known that this is a region of the 
atmosphere in which in general the mesoscale models are 
less sensible than near the ground. It is therefore a further 
indication that the Meso-Nh model is well placed to forecast 
the turbulence evolution at these time scales. 

Once the tendency in overestimating the strength of the tur- 
bulence in the surface layer will be solved (forthcoming pa- 
per) we plan to run the Meso-Nh model in other regions 
of the internal Antarctic Plateau to identify the best loca- 
tions for astronomical observations i.e. the places with the 
best turbulence characteristics from an astronomical point 
of view. 
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